function Hv = H( a, b )
%H Возвращает вектор фокусировки для нашей антенной системы
%   Входные аргументы - фзимут и угол места

[kx, ky, kz] = sph2cart(b, a, 1);
k = [kx;ky;kz];

r{1} = [1/4, -1/4, 0];
r{2} = [-1/4, -1/4, 0];
r{3} = [-1/4, 1/4, 0];
r{4} = [1/4, 1/4, 0];
NAnt = length(r);

phi = nan(NAnt, 1);
for i =1:NAnt
    phi(i) = 2*pi * r{i}*k;
end

Hv = exp(1i*phi);
end

